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Abstract 

We investigate a model of large extra dimensions where the internal space has the geometry of a 
hyperbolic disc. Compared with the ADD model, this model provides a more satisfactory solution 
to the hierarchy problem between the electroweak scale and the Planck scale, and it also avoids 
constraints from astrophysics. In general, a novel feature of this model is that the physical results 
depend on the position of the brane in the internal space, and in particular, the signal almost 
disappears completely if the brane is positioned at the center of the disc. Since there is no known 
analytic form of the Kaluza-Klein spectrum for our choice of geometry, we obtain a spectrum based 
on a combination of approximations and numerical computations. We study the possible signatures 
of our model for hadron colliders, especially the LHC, where the most important processes are the 
production of a graviton together with a hadronic jet or a photon. We find that the signals are 
similar to those of the ADD model, regarding both qualitative behavior and strength. For the case 
of hadronic jet production, it is possible to obtain relatively strong signals, while for the case of 
photon production, this is much more difficult. 
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I. INTRODUCTION 



The Large Hadron Collider (LHC) at CERN near Geneva, Switzerland is about to be- 
come operative. The searches at the LHC for new physics beyond the Standard Model (SM) 
will mainly include the potential discoveries of the Higgs boson, supersymmetry, and extra 
dimensions. In this paper, we will be interested in the third issue, i.e., extra dimensions. In- 
deed, an observation of extra dimensions would be truly revolutionary and would completely 
change our view of the Universe. 

The idea that spacetime could have more than four dimensions was first proposed by 
Theodore Kaluza [l] and Oskar Klein at the beginning of the twentieth century. One of the 
most interesting features of extra dimensions is that they are not ruled out by experiments, 
provided only that they are compact and small enough to have avoided detection so far. In 
the scenario known as large extra dimensions, they could even be macroscopically large. 

Large extra dimensions were first proposed by Arkani-Hamed, Dvali, and Dimopoulos 
(ADD) in 1998, their model being known as the ADD model 3|, |4j. The novel feature of this 
model is the assumption that the SM fields are confined to a so-called brane, which is a four- 
dimensional manifold residing in the full bulk spacetime. This brane is to be identified with 
ordinary four-dimensional spacetime. Since the SM fields are not allowed to probe the extra 
dimensions, experimental constraints on their size are avoided to a large extent. Gravity, 
on the other hand, carries no SM charges and is allowed to probe the extra dimensions. 
In principle, the assumption that gravity lives in a higher-dimensional spacetime leads to 
sizable deviations from Newton's inverse-square law at short distances. However, because 
of the weakness of the gravitational force relative to the SM forces, Newton's law has only 
been tested down to distances of the order of micrometers, and hence, the experimental 
constraints are still quite weak. 

One of the main motivations for the ADD model is that it provides a solution to the so- 
called hierarchy problem between the electroweak scale M ew ~ 100 GeV and the (reduced) 
Planck scale Mpi ~ 10 18 GeV. Theoretically, the bare Higgs mass is expected to receive 
higher-order quantum corrections of the order of Mp\. This would mean that extreme fine- 
tuning of the parameters would be needed in order for the electroweak scale to be as low as 
100 GeV. In fact, the ADD model provides a very elegant solution to this problem. Since 
gravity really propagates in more than four spacetime dimensions, the Planck scale that 



2 



we observe through gravitational measurements is an effective scale, valid only for energies 
lower than the inverse of the radius of the internal space. The Planck scale is related to the 
true fundamental energy scale for gravity through the volume of the internal space. If this 
volume is large enough, then the fundamental scale for gravity could actually be as low as the 
electroweak scale. However, there is a problem related to this solution in the ADD model. 
While the problem of the hierarchy between the electroweak scale and the fundamental scale 
for gravity is solved, there is a new large hierarchy between the electroweak scale and the 
inverse of the radius of the internal space. Thus, the hierarchy problem is only reformulated 
as the question of why the radius of the internal space is so large compared to the electroweak 
scale. 

However, in the ADD model, the internal space is assumed to be flat and compactified on 
a torus. Thus, one possible solution to the problem of the hierarchy between the electroweak 
scale and the radius of the internal space is to drop this assumption and instead consider 
a different geometry. Therefore, it has been argued in Ref. [5] that a compact hyperbolic 
internal space in particular is a better alternative than the flat geometry of the ADD model. 
Note that, in some sense, the hyperbolic model is a generalization of the ADD model. 

In addition, it should be mentioned that there are other models of extra dimensions that 
include branes. One of the most important models is the so-called Randall-Sundrum (RS) 

n q 

model [a, 13|) in which two branes are introduced and the SM fields are confined to one of 
these branes only. Nevertheless, we will not consider such models further. 

In this paper, we investigate large extra dimensions with the internal space being a two- 
dimensional hyperbolic disc. Especially, we study two plausible signals, i.e., the reactions 
p + p — > jet + G and p + p — > 7 + G, where G denotes a Kaluza-Klein (KK) mode of the 
graviton, that could be measured at the LHC using missing-energy techniques. It should be 
noted that the hyperbolic disc model, like the ADD model, is only an effective theory, which 
means that it is a non-renormalizable low-energy approximation of a more fundamental 
theory that is called the ultraviolet (UV) completion of the effective theory. 
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renomenology of the ADD model has been extensively investigated in the literature 
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1 1 51 ] - In particular, signals of the ADD model that are relevant for the 
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LHC have been studied in Ref. 



. In Ref. 



16] , a model with a spherical internal space has 



been examined. In addition, a model of RS type, which is similar to ours, was considered in 
Ref. [ill , in the setting of discretized extra dimensions. Finally, hyperbolic extra dimensions 
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could have interesting implications in cosmology, which have been studied in Refs. [18|, ll9| , 
though we will not discuss this issue further in this paper. 

This paper is organized as follows. In Sec. [TTJ we present the hyperbolic disc model 
for large extra dimensions and obtain an approximate form for the KK spectrum of the 
graviton in this model. Then, in Sec. IHH we analyze the interactions between the graviton 
and the SM fields that are relevant for the plausible signals of the model at the LHC Next, 
in Sec. IIVt we give our numerical results for the cross sections of the signals discussed in 
Sec. IIHI Finally, in Sec. [V] we summarize our results and present our conclusions. 



II. THE HYPERBOLIC DISC MODEL 
A. Hyperbolic extra dimensions 

The model that we consider is similar to the ADD model, with the only exception that 
the internal space is a two-dimensional hyperbolic disc, which is denoted H 2 . Hence, the 
geometry of the higher-dimensional spacetime is a product M 4 x H 2 , where M 4 denotes 
four- dimensional Minkowski space. The SM fields are assumed to be confined to a four- 
dimensional brane, while gravity alone probes the extra dimensions. The metric for the 
six-dimensional spacetime is 

(9mn) = diagfl, -1,-1, -1, -1, -v~ 2 sinh 2 (wr)], (1) 

where r G [0, L] and ip G [0, 2tt) are polar coordinates and v is the curvature of the disc. 
The coordinate system is such that r is the physical radial distance between the origin and 
a point (r,<p). We follow the convention that indices in the full spacetime are written as 
upper-case Roman letters, M = 0, 1, 2, 3, 5, 6, Minkowski indices are written as lower-case 
Greek letters, fi = 0, 1, 2, 3, and indices in the internal space are written as lower-case Roman 
letters, i = 5,6. Also, x denotes the coordinates in M 4 or in the full higher-dimensional 
spacetime and y the coordinates in H 2 . Note that |det(<7Miv)| = |det(<7y)|, which means 
that there is no ambiguity in using the symbol \g\ for both of these quantities. The number 
of extra dimensions is denoted by d. 

The most common way to hide the extra dimensions is through compactification of the 
internal space as a quotient space H 2 /T, where T is a discrete subgroup of the isometry 
group of the internal space. In this paper, we consider instead an internal space with an 
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explicit boundary. The main motivation for this choice is computational simplicity. For 
an internal space of hyperbolic geometry, it is not possible to solve for the KK spectrum 
analytically. Instead, numerical calculations are needed, and these are much simpler in a 
space with a boundary than in a quotient space. We do not attempt to describe the origin of 
the boundary, but simply to investigate its possible implications. An important consequence 
of this choice of geometry is that, in contrast to the ADD model, the physical results depend 
on the position of the brane in the internal space. 

The energy-momentum tensor corresponding to the metric ([1]) is {Tmn) — 
diag(— v 2 , v 2 , v 2 , v 2 , 0, 0). As mentioned, we do not attempt to justify this solution of Ein- 



stein's equations. For a deeper discussion, see Refs. [171. |20|. 

The most important advantage of a hyperbolic space is that it offers the possibility of a 
more satisfactory solution to the hierarchy problem than the ADD model does, as is described 
below. Another important advantage is that astrophysical constraints on the lower bound 
on the fundamental mass scale, which are particularly important in the two-dimensional 
ADD model 2lj, can be avoided to a large extent. Thus, the model can allow for a low 
value of this mass scale even in the case of two extra dimensions only. The constraints on 
the parameter space are described in more detail in Sec. Ill Dl 

Since gravity is the only field probing the internal space, it plays an important part in any 
phenomenological studies of the model. By assumption, it is governed by the six- dimensional 
Einstein-Hilbert action 

where R is the Ricci scalar, A is a cosmological constant, and the mass scale M* is introduced 
in order to make the action dimensionless. In the same way as in the ADD model, M* 
replaces the Planck scale Mpj as the fundamental mass scale for gravity. The two scales are 
related through the equation Mpj = VM* 21], where V is the volume of the internal space. 
However, note that the definition of the fundamental mass scale differs between authors. In 
order for the model to provide a solution to the hierarchy problem, we demand that M* 
is of the order of 1 TeV. In the ADD model, where Vd = (2nL) d , this gives the radius 
L ~ io 31 / d TeV -1 , which is unnaturally large in comparison to if d is not very large. 
Hence, the hierarchy problem is not really solved, but simply rephrased as the question of 
why the product M*L is large. In our model, on the other hand, the volume of the internal 
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space is 

V = J dV = J^dr ^^v^j=^sinh 2 {^j ■ (3) 

For large vL, the volume increases exponentially as a function of the radius. This has the 
consequence that the relation between the mass scale and the volume of the internal space 
can be satisfied for M* ~ 1 TeV without generating a large hierarchy between M* and L, if 
v is suitably adjusted. This fact is our main motivation for the hyperbolic geometry of the 
internal space. In Fig. [Tj, the product M*L is plotted as a function of v. The result is nearly 
independent of M* in the range that we are interested in. For v = M*, we obtain the lowest 
possible value M*L ~ 100, while for smaller v, the value of the product is significantly larger. 
Note that our effective model is supposed to be valid only up to energies of the order of M*, 
and hence, we do not consider values of v larger than this scale. Thus, the best possible 
solution to the hierarchy problem in our model is obtained when the curvature v is of the 
same order of magnitude as the fundamental mass scale. In this case, there is also no new 
hierarchy problem involving v. 




FIG. 1: The product M*L function of v/M*. 



B. Kaluza Klein decomposition of the graviton 

In order to investigate the phenomenology of our model, we now use the ordinary proce- 
dure of KK decomposition to reformulate it as an equivalent four-dimensional field theory. 
In this picture, the graviton field propagating in the full spacetime is represented by an 
infinite KK tower of particles with different masses, called the KK modes of the graviton. 
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The set of masses is known as the KK spectrum, and each mass corresponds to a quantum 
of momentum in the internal space. 

In the ADD model with d extra dimensions, the graviton living in the full spacetime is 
expanded in terms of harmonic functions exp (in • y/R) with corresponding masses m n = 
\n\/R. Here, n is a (^-dimensional vector with integer entries. In our model, the geometry is 
more complicated, and this has the result that it is not possible to obtain the KK spectrum 
analytically. Thus, in this section, we derive approximate expressions for the KK modes and 
masses in our model. 

The starting point is the equations of motion for the graviton. From the action (T5]), it 
follows that the free equations of motion are the free Einstein equations, i.e., Rmn — 0, 
where Rmn is the Ricci tensor. The dynamics of the graviton is studied by expanding 
these equations to first order in a perturbation Iimn about the background metric ([1]). The 
perturbation h MN is interpreted as the massless spin-2 graviton field in the six-dimensional 
spacetime. We are only interested in the dynamics of the four-dimensional part h^ v of 
the perturbation, and hence, we make the simplifying approximation of setting all other 
components to zero, which has often been done in the literature fltl 17 . 22]. Therefore, the 
resulting metric is 

( , (vr> + K u {x,y)/M* \ 
g M N{x,y) - , (4) 

V 9iM ) 

where rj^ is the Minkowski metric and the factor M~ 2 ensures that the graviton has the 
correct dimension, i.e., (mass) 2 . Under the assumption that the unperturbed Einstein equa- 
tions are satisfied, the detailed derivation of the linearized equations of motion has been 
performed in Ref. JliJ. The resulting equations are 

A LB V = 0, (5) 

where Alb = V m Vm is the Laplace-Beltrami (LB) operator, which is the generalization of 
the Laplace operator to curved spaces. In terms of a coordinate system {x M }, 

1 „ / rr^MNi 



A LB ^ = ^=d M ( ^/\g\g MN d N ij ) . (6) 

g\ 



Because of the factorizable geometry, the Laplace-Beltrami operator can be written as 
Alb = D + A#2, where □ = d^d^ is the d'Alembert operator in four- dimensional Minkowski 
space and A#2 is the Laplace-Beltrami operator in the two-dimensional hyperbolic space. 



In order to reformulate the theory without explicit reference to the extra dimensions, h^ v is 
expanded in terms of the eigenfunctions of A#2, i.e., 

h^(x,y) = ^2K,iw(x)if> n (y). (7) 

n 

The coefficient functions h n>fU ,(x) are the KK modes of the graviton, satisfying the equations 

(□ + m 2 n ) h n ^ u = 0, (8) 

where m 2 n is the eigenvalue corresponding to the eigenfunction ip n . Here, n denotes any 
general set of Kaluza-Klein indices. From Eq. (jSJ), it follows that m n has the interpretation 
of the mass of the KK mode h n ^ u . 

Although no analytic form for the KK spectrum is known, several important results are 
generally true for the spectrum of the Laplace-Beltrami operator on a Riemannian manifold 
M (possibly with a boundary). It is assumed that the closure of M is connected and 
compact, which is true for a hyperbolic disc. The eigenfunctions belong to the Hilbert space 
L 2 (M) of square- integrable functions on M, with the inner product given by 

(/,#> = J d d yf(yyg(y), (9) 

where the star denotes complex conjugation. Then, if the boundary conditions fall into one 
of the four categories a) Dirichlet conditions, b) Neumann conditions, c) mixed Dirichlet 
and Neumann conditions, or d) periodic conditions, the following results hold (23! ]: 

1. The set of eigenvalues consists of a sequence, = Ai < A2 < • . . 00 in the case of 
Neumann or periodic conditions, or < Ai < A 2 < . • . 00 in the other cases, and each 
associated eigenspace is finite-dimensional. 

2. Eigenspaces, belonging to distinct eigenvalues, are orthogonal in L 2 (M), which is the 
direct sum of all eigenspaces. 

Thus, when the internal space is compact, the spectrum of the Laplace-Beltrami operator 
is countable, and it is possible to make a KK expansion of the type (j7]). 

Now, we need boundary conditions at r = and r = L. At r = 0, we only have to demand 
that the solutions are finite. At r = L, the possible alternatives are in principle Dirichlet 
or Neumann conditions, or a combination of both. As mentioned above, the eigenvalues are 
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interpreted as the squared masses of the corresponding KK modes, i.e., X n = m 2 . According 
to[TJ, there is no massless mode in the spectrum for Dirichlet or mixed conditions, while for 
Neumann conditions there is such a mode. Hence, in order to obtain the correct low-energy 
behavior, we impose Neumann conditions at r — L. 

We also make use of the following result regarding the asymptotic distribution of the 
eigenvalues, known as Weyl's asymptotic formula: If N(X) denotes the number of eigenvalues 
in the interval [0, A], counted with multiplicity, then 24| 

N(X) u d V d 



lim 



(10) 



A-,00 \d/z (27r) d ' 

where d is the dimensionality of the manifold M, Ud is the area of the unit disc in M d , and 



Vd is the volume of M. Taking M to be the hyperbolic disc, and using Eq. 
the result 



lim 



N(m 2 ) sinh 2 (t>L/2) 



m-+oo m 2 v 2 



we obtain 



(11) 



We now find the general solution of the eigenvalue equation A^ip = m 2 if). These eigen- 
functions contain important information on the coupling of the graviton KK modes to SM 
fields, and they are also the starting point for our numerical investigations of the KK spec- 
trum. For the case of hyperbolic geometry, the Laplace-Beltrami operator is given by 



1 







sinh(tr) dr 



sinh(vr 



dip 
dr 



d 2 ip 



sinh (vr) dip 2 



(12) 



Introducing the dimensionless parameter r = vr, we obtain the eigenvalue equation 



1 



d 



sinh(r) dr 



sinh(r 



dip 



1 



d 2 iP 



sinh (t) dip 2 



k 2 iP, 



(13) 



where k 2 = m 2 /v 2 . In order to solve this equation, we first expand ip in the angular direction 



Hip 



(14) 



and find the radial equation 



sinh(r)— — 
dr 



sinh(r) dr 

Introducing x = cosh(r), we have the equation 



sinh (t) 



-To = k 2 T f . 



(15) 



d_ 

dx 



[l-x 1 



dx 



+ 



u(v+l) - 



Tt = 0, 



x' 



(16) 
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where u(u+l) = —k 2 . This is Legendre's associated equation. Its solutions are the associated 



Legendre functions of the first and second kind [25(. Note that these are not the same as 
the functions encountered e.g. in spherical harmonics, since the domain here is [1, oo) rather 
than [—1, 1]. In particular, in this case, v is not restricted to integer values in the interval 
[—£,£}. The associated Legendre functions of the second kind are divergent at x = 1, and 
can therefore be discarded as non-physical solutions. The associated Legendre functions of 
the first kind are the physically acceptable solutions. In what follows, they will be referred 
to simply as the Legendre functions and they will be denoted by P„. They can be expressed, 
up to normalization, as 

Pt(x) = W2 F (-„, u+l,l + \l\, , (17) 

where F is the hypergeometric function. It is common to introduce the parametrization 
v = — 1/2 + ip, since for real p, P_ x i 2+ip is real. Note that, since v = — l/2±ip both give the 
same value for v(v + l), we only need to consider values of p in e.g. the half-plane Re(p) > 0. 
In terms of the parameter p, the eigenvalues are 

m% = v 2 (± + f?). (18) 

Finally, the radial eigenfunctions are 

"1 . 1 



T» = P\ +ip [cosh(vr)} = tanh"« (^) F - - ip, - + ip, 1 + \£\, - sinh 2 



(19) 



An important consequence of this result is that, since F(a, b, c, 0) = 1, T p i(0) = for I ^ 0. 
Hence, most of the eigenfunctions are equal to zero at the origin. As is demonstrated in 
Sec. IIII Al the couplings of the KK modes to SM fields are proportional to the modulus 
squared of the eigenfunctions, evaluated at the position of the brane. Thus, for the most 
symmetric location of the brane, at r = 0, only the I = modes couple to SM fields. In this 
case, it is not possible to probe the extra dimensions with the methods that are investigated 
in this paper, as is discussed in Sec. HVl 

If the parameter p is restricted to real values, then the relation (1181) implies that the 
spectrum is restricted to the interval [v/2, oo). However, there is a priori nothing to prevent 
p from being complex. The spectrum of the Laplace-Beltrami operator is real and non- 
negative, and m = is in one-to-one correspondence with the constant eigenf unction, but 
there could possibly exist eigenvalues in the interval (0, v/2), corresponding to values of p 
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in the imaginary interval (0, i/2). We have numerically investigated the zeros of dTpi/dr for 
p in this interval and found none. However, we have not been able to prove this result, and 
such a proof would, of course, be of interest. Nevertheless, we assume in the remainder of 
this paper that the KK spectrum lies in the interval [v/2, oo). Thus, under this assumption, 
there is a mass gap between zero and the mass m\ = v/2 of the first KK mode. The 
significance of this result is discussed in Sec. Ill Dt 

As eigenfunctions, the normalization of the functions T p t is not determined. It is decided 
by the normalization of the Lagrangian kinetic terms of the individual KK modes. As the 
higher- dimensional kinetic terms involving derivatives with respect to r and tp become mass 
terms in the four- dimensional picture, those terms are irrelevant for the following discussion. 
Since we will later use results for the ADD model from Ref. jg|, we follow their convention, 
where the relevant kinetic terms are of the forms 

4 kin) = \d x h^-AK^. (20) 

Furthermore, since, in the ADD model, h pu (x,y) = J2n hn,pu(x) exp (in • y/R) and is 
real, we must have h^ n = h£ u * , which means that 

4 km) = \d x hrd x hn,,u. (21) 

In our model, we have 



S (khl) = f V\9\d 6 x^d x h^d x h, u 



J d 4 x j >/W 2 v\& (E^^J ^ ^EV',pV 

p,i p',£' J ~ ■> 

EE / ^^^hp^d^p't'^pp'S-zA^pA? 

p,£ p',£' J 

E / d A x\M\ 2 ^d x h^_ e d x h pe ^ 

E / d 4 s||Vyll a ^r^W> (22) 



2 

P ,t " 

where in the last equality we have used the fact that ip p -t = which follows from the 
results (T p -i)* = T p ^g = T P £ and exp (i£(p)* = exp (— i£<p). Since h pu is real, this result 
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implies that h^_ t = . The rather than a &w, in the fourth line comes from the 

fact that there is no complex conjugation on i/j p e in the third line. Thus, the normalization of 
the individual kinetic terms is the same as in Eq. (1211) . if we set ||Vvl| 2 = ^> which determines 
the overall normalization of the eigenfunctions. 



C. Approximate eigenfunctions 

In principle, all the information about the eigenfunctions is given by Eq. (fl9l) . However, 
in order to better understand their behavior, it is useful to consider a certain approximation 



of them, which has been adapted from a similar case in Ref. 261, and is based on the 

n 

Wentzel-Kramers-Brillouin (WKB) approximation [27J, familiar from quantum mechanics. 
The approximate expressions for the eigenfunctions given by this approximation also have 
the advantage that their numerical evaluation requires significantly less computer power 
than the exact expressions, which has been important for the calculations in Sec. IIVI In 
order to find these approximate expressions, we introduce the auxiliary functions u p i(t) = 



a/ sinh(r)T p ^(r). In terms of u p e, Eq. ( fl5]) becomes 

d 2 u p£ I 2 - 1/4 2 
~ ^ + Upi = P Upl - (23) 

This equation has the form of a one-dimensional Schrodinger equation with energy E = p 2 
and potential V(r) = (£ 2 — 1/4)/ sinh 2 (r). Its solutions have differing qualitative behavior 
depending on the relative magnitudes of E and V(t). For E > V(t), they are oscillatory, 
while for E < V(r), there is one increasing and one decreasing solution. The turning point To 



between the two regions, given by the equation E = V(r ), is r = arsinh \/ (£ 2 — 1/4)/ p 2 
For r < 1, the solutions are approximately u p e(r) = sinrr^'tY). The decreasing solutions 
diverge at the origin, and correspond to the Legendre functions of the second kind, while 
the increasing ones correspond to the Legendre functions of the first kind. 
Now, for r > r , the WKB approximation gives the solutions 

. s sin[6(r)l . . 

VW = " 1 J ll/4 > r>r , (24) 

2 l 2 -l/4 1 1 



where 




0(r)= / dr\lp 2 --^^pr + ^. (25) 
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Here, we have also used the fact that the eigenfunctions are real. For r < tq, T p i is small and 
is approximated as zero, with the understanding that there are no zeros of the derivative in 
this region. In order for the eigenfunctions to be continuous at r , we set the phase <p = pr . 
Thus, the approximate expression for T p i that we use is given by 

{sm[p(r-T )] r > r 

[pW(t)-(*3-i/4)] 1/4 ' - (26) 
0, r < Tq. 

The condition for the WKB approximation to hold is 

dA 



1 
2^ 



dr 



« 1, (27) 



where 



2n 27rsinh(r) 
(r) = [E — V(r)] 1 / 2 = [p2 S inh 2 (r) - (£ 2 - ( ^ 

This function blows up at r = tq, where the WKB approximation is generally not valid. 

A sample of the approximate functions, as well as the corresponding exact functions, 
are presented in Fig. [2j We have plotted the squared absolute values of these functions, 
since these are the quantities that enter in the physical results. For small vr, where the 
exact functions have not yet started to increase appreciably, the approximations agree with 
the exact functions to good accuracy. As expected, the agreement is worse in the region 
around t = tq, where the approximate functions start to oscillate. From the first minimum 
in the oscillating region and on, the approximation is once again very accurate, and the 
agreement increases with increasing vr. Also, our numerical investigations indicate that the 
approximation becomes better as the parameters p and t are increased. 



D. Constraints on the parameter space 

In the ADD model, there is only a single free parameter, which can be taken to be the 
fundamental mass scale M*. In our model, on the other hand, the curvature v of the internal 
space and the position of the brane in the radial direction tj, enter as two additional free 
parameters. The mass scale M* is bounded from below to M* > 1 TeV by the fact that 
no signs of quantum gravitational effects have been found in experiments up to this scale 



28( . On the other hand, M* cannot be much larger than 1 TeV if the model is to provide a 



solution to the hierarchy problem. The fact that our model is an effective theory, valid only 
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FIG. 2: The squared absolute values of the approximate eigenfunctions (j26[) compared to the exact 
eigenfunctions (|19p as functions of vr. The solid curves are the exact functions and the dashed 
curves are the approximate functions. In the left panel, the functions are plotted for the parameter 
values p = 5 and £ = 2, whereas in the right panel, they are plotted for the parameter values p = 2 
and i = 50. 

up to energy scales of the order of M*, means that we should not consider values of v larger 
than this scale. 

Further constraints can be found by demanding that the model should not be in conflict 
with other well-established physical phenomena. In the context of the ADD model, a num- 
ber of such constraints have been analyzed in Ref. 
come from astrophysics and cosmology. If sufficiently light, the lightest KK mode could be 
produced in large numbers in high-temperature systems, such as supernovae, and carry away 
large amounts of energy. This could potentially alter the evolution of the system in a non- 
acceptable way. For the ADD model, this places important constraints on M* In our 
model, the mass of the lightest KK mode is bounded from below by mi = v/2. If v is chosen 
so that this mass is larger than the temperature of a supernova, which for SN1987A is about 
50 MeV, then the constraints are completely avoided. This is achieved for v > 100 MeV. 
Note that this bound has not been optimized, but merely gives an order-of-magnitude es- 
timate. Of course, the bound also depends on M*. However, a more detailed analysis of 
the exact constraints on the full parameter space is beyond the scope of this paper. As has 
been mentioned earlier, we only consider values of v of the order of M*, in order to obtain 
a satisfactory solution to the hierarchy problem. Since M* is of the order of 1 TeV, v is far 
larger than 100 MeV. In the same way, we consider r& to be unrestricted to lie anywhere in 
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the range [0, r max ], regardless of M* and v. 



E. Numerical analysis of the Kaluza Klein spectrum 

Since it is not possible to obtain the KK spectrum of the graviton analytically, we have 
used a combination of numerical calculations and Weyl's asymptotic formula, as well as the 
form of the approximate solutions (126]) obtained by using the WKB approximation. We 
assume that there are no zeros of the derivative of T p ^{r) for r < r . Thus, for a given value 
of £, this means that any allowed value of p has to be such as to fulfill the relation 



where r max = vL. 

Performing the numerical calculations, we have found that the spectrum increases loga- 
rithmically for small m, i.e., m n ~ log(n). Comparing it with the result of Weyl's asymptotic 
formula, which increases as m n ~ ^/n, it is in fact nearly constant at m w mi = v/2. This 
is expected, since the spectrum from Weyl's formula starts out at m = 0, while the true 
spectrum starts out at the non-zero value mi. Hence, the true spectrum has to increase 
slower than the approximate formula in order for the two results to converge for large n. 
We have not been able to solve numerically for the spectrum up to values where the two 
results converge. Instead, we have resorted to solving numerically only for a manageable 
number, (9(10 3 ), of eigenvalues and extrapolating these results up to a point where Weyl's 
formula is supposed to hold. This point is taken as the intersection between the extrapolated 
results from the numerical solution and the result from Weyl's formula. Since the spectrum 
is nearly constant at m = mi in the lower regime, our final result is 



where Q(x) is the Heaviside step function. 

III. INTERACTIONS BETWEEN THE GRAVITON AND THE SM FIELDS 
A. The interaction Lagrangian 

In order to study interactions between the KK modes of the graviton and the SM fields, 
we need the interaction terms in the action. In the general case when tiMi 7^ 0, the higher- 




(29) 




(30) 
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dimensional coupling between gravity and the SM fields to first order in Hmn/M^ is given 
by 

^ (mt) = J^J dAx J d2 y TMN ^y) h MN(x,y)- (si) 

Because of the confinement of the SM fields to the brane at y — yb, the energy-momentum 
tensor is T MN (x,y) = 6jf S^T^ U (x)S^ (y — yb) (no summation), where S^ 2 '(y) is the Dirac 
delta function and T^ u {x) is the ordinary SM energy- momentum tensor. Inserting this 
expression into Eq. (131]) and using the KK expansion ((Tj) yields the individual interaction 
terms 

tft* = ^J d 4 *T^(*)^V,^)^), (32) 

where we have explicitly displayed the normalization of • Using the relation Mp } = VM^, 
these interaction terms can be rewritten as 

S ''°" = JftJ d4 * T ^MWM, (33) 



(34) 



where 



which are dimensionless numbers characterizing the coupling strengths. In the ADD model, 
the corresponding constants c n = exp(in-y/R) are unimodular, and hence, they do not 
affect any physical results. In our model, the numbers c p e are generally not unimodular, 
which has the results that different KK modes have different coupling strengths and that 
these coupling strengths depend on the position of the brane in the internal space. However, 
note that the angular parts of the eigenf unctions, exp (i£(p), are still unimodular, and thus, 
Cpi only depends on the radial position of the brane and not on the angular position. From 
the above discussion, it follows that the constants c p i provide a parametrization of the 
difference between the ADD model and the hyperbolic disc model in the coupling strengths 
of the individual KK modes. 



B. Graviton production cross sections 

Since the interaction terms in the action (133]) differ from those in the ADD model only 
by the constant factors c p e, the Feynman rules for the hyperbolic disc model are the same as 
those for the ADD model, except that the vertex factors involving KK modes of the graviton 
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are multiplied by these same factors. The Feynman rules for the ADD model are given in 
Ref. 8|, and since our normalization convention agrees with theirs, the results can be used 
to obtain the Feynman rules for the hyperbolic disc model. Note, though, that their results 
are expressed in terms of the quantity Mq, which is related to M* as = \/2irM*. 

Since the amplitudes for processes involving KK modes of the graviton are suppressed 
by the Planck scale, the cross sections for production of a single KK mode are extremely 
low. However, at sufficiently high energies, a large number iV(m^ ax ) = N(E 2 m ) « 
v ~ 2 sinh 2 (vL/2) E 2 m of KK modes are kinematically available. The cross sections for pro- 
duction of any available KK mode are suppressed only by powers of the higher-dimensional 
gravitational mass scale M*, which could be significantly smaller than the ordinary Planck 
scale. Once produced, a KK mode is extremely weakly interacting and consequently it 
appears as missing energy in detectors. 

The cross section for production of any KK mode is obtained by summing over the 
kinematically available individual cross sections, i.e., 

da \ ~\ da m 

where s and t are the usual Mandelstam variables. The standard way to treat this sum is 
to observe that the mass splittings are small in relation to any other relevant mass scale. In 
the ADD model, the mass splittings are Am ~ R~ l , while in our model, the mass splittings 
are Am ~ wsinh -1 (vL/2), which is small for the regions of the parameter space that we 
consider. Thus, the sum can be approximated to good accuracy by an integral, i.e., 

— — « / dmn(m)— ^, (36) 
dt J K ' dt y J 

where n(m) is the density of states for the KK modes. Using the approximate expression 
for the spectrum given in Eq. (1501) . the density of states is 



Tt ~ ^ "dP (35) 



dN(m 2 ) r/ N sinh 2 (wL/2) 2 n , . 2 sinh 2 (uL/2) 

— — L = S(m - mi) v 1 ' m 2 + Q(m - m{) \ — 

dm v 2 v l 



There is a technical problem related to the evaluation of the integral (|36l) . While the 
differential cross section da m /dt and the density of states n(m) are given as functions of 
m, the constants c p e are only available as functions of the KK indices p and £, and we have 
no analytic relation between these quantities. Hence, in order to evaluate the integral, we 
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need to numerically translate the constants into functions of m. This can be done by av- 
eraging c p e over the allowed values of I for each value of p, and using the relation (|18p to 
express p in terms of m. From Eq. f!29p . it follows that, for fixed p, i is restricted to the 
interval [— £ max , £ m ax\, where £ max = \J p 2 sinh 2 (r max ) + 1/4. Now, in order to perform the 
averaging, an estimate of the density of eigenvalues for given p and I is needed. The ap- 
proximate eigenfunctions (|26|) that we employ consist of a decaying factor and an oscillating 
factor sin[p(r — r )]. Since the eigenvalues are determined by the positions of the zeros of 
dT p i/dr\ T=Tb as a function of p, we consider the eigenfunctions evaluated at the position of 
the brane, i.e., at r = r^. In addition, since tq = arsinh 



is a slowly varying 



function of p and £, we may locally consider the oscillating factor, seen as a function of p, to 
have a well-defined wave number equal to r b — r . Thus, the density of zeros of T p i(r b ) as a 
function of p is approximately proportional to r& — tq. We are interested in the zeros of the 
derivative dT p i/dr\ T=Tbl and we assume that there is exactly one such zero between each pair 
of zeros of T p £ (r^). Hence, for given p and £, the density of eigenvalues is also proportional 
to T b — t . Using this result, the averaging has been performed. A sample of the resulting 
functions is shown in Fig. [3j Note that in the resulting coefficients c m , the rapid oscillations 
of the eigenfunctions have been washed out, leaving smooth functions. 
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FIG. 3: The squared absolute values of the averaged coefficients \c m \ 2 as functions of m. The black 
curves correspond to the parameter value r& = r max and the gray ones correspond to the value 
T~b = Tmax/2- As a reference, the corresponding trivial results for the ADD model are also shown as 
dotted lines. In the left panel, the results are plotted for M* = 1 TeV, whereas in the right panel, 
they are plotted for M* = 2 TeV. 
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Finally, we need to take into account that the colliding particles are protons, while the 
cross sections are given on the more fundamental quark level. Cross sections for such pro- 
cesses are calculated using the parton model. The total cross section for a high-energy 
hadron-hadron collision can be written as a convolution of two parton distribution functions 
with a hard-scattering parton- level cross section a 30] , 

a A +B-,x{s) = ^ / dxi / dx 2 f a /A{xi,s)fb/B{x2,s)cra+b^x{s), (38) 

0,6 J ° J ° 

where f a /A is the parton distribution function for the parton a in the hadron A, s is the 
squared center-of-mass energy in the hadron-hadron system, x\ and x 2 are the momentum 
fractions of the hadrons carried by the respective partons, and s = X\X 2 s is the effective 
center-of-mass energy squared in the parton-parton system. In principle, the sum is to 
be taken over all parton species, i.e., quarks, anti-quarks, and gluons, although the heavy 
quarks, i.e., charm, bottom, and top, are usually neglected. In this paper, we use the 



CTEQ6M 



3JJ parton distribution functions. 



C. LHC graviton production reactions 

We consider the reactions p + p — > jet + G and p + p — > 7 + G, where G denotes a KK 
mode of the graviton. For the ADD model, the individual differential cross sections for these 
reactions are given in Ref. On the parton level, the reaction p + p — > jet + G consists 
of the three subprocesses q + q — > g + G, q + g — > q + G, and g + g — > g + G, while the 
reaction p + p — > 7 + G consists of the single subprocess q + q — > 7 + G. Since each of these 
subprocesses includes a single vertex involving a graviton, the cross sections for our model 
are obtained by multiplying the results for the ADD model by the constant factors \c p t\ 2 . 

As mentioned above, the produced gravitons are very weakly interacting, and hence, they 
appear as missing energy in detectors. Thus, the observed reactions are p + p — > jet + f and 
p + p — > 7 + respectively, where ft denotes the missing energy. 



IV. NUMERICAL ANALYSIS OF SIGNALS 



In this section, we present the predictions of our model. For both cases of jet and photon 
production, we have calculated differential cross sections with respect to cos(#), where 9 is 
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the angle between the proton beam and the outgoing jet/photon, as well as with respect to 
PrjetAp which is the momentum of the jet or photon perpendicular to the beam. All of our 
results are presented for yfs = 14 TeV. 

Since we consider high-energy processes, the partons are approximated as being massless. 
As there are massless particles in the final states in both of the processes that we study, 
the expressions for the cross sections suffer from collinear divergences in the limit of zero 
transverse momentum of these particles. In order to avoid these singularities, we impose 
a lower cut-off on the transverse momentum, pr,7/jet > p^ m - This cut-off also serves to 
increase the signal-to-background ratio. Because of the finite size of the detector, there 
is an upper cut-off on the longitudinal rapidity (or pseudorapidity) i] = artanh [cos(#)] of 
these outgoing particles, i.e., |n 7 /j e t| < i] mauX . We have used the value ?7 max = 2.5 for all 



measurements considered 



32]. 



It is important to take into consideration the fact that the theory is an effective one 
only, which is supposed to break down at large energies, above some cut-off scale. Without 
knowledge of the UV completion of the effective theory, it is not possible to determine this 
cut-off scale exactly. We follow Ref. {s] and trust our results only up to the mass scale 
Mu = \/2txM*. We also follow their method of analyzing the validity of the results, i.e., 
by computing cross sections that are set to zero for s > M|>, and compare these to the 
naive results. In regions where the results agree, almost all of the contributions to the 
cross sections come from subprocesses with an effective center-of-mass energy lower than 
the fundamental mass scale, and hence, these results can be trusted. These regions depend 
on the chosen set of parameters, and differ between jet and photon reactions. In general, 
the results become better for higher M*, but at the same time the cross sections decrease. 
Hence, we need to make a trade-off between these two competing effects. 

As discussed in Sec. [Til we only consider values of v of the order of M*. However, for 
v = M*, it is difficult to to obtain valid results for the effective model. Hence, for all the 
results presented, we have have set v = M*/2, in which case it is possible to obtain valid 
results. Complementary to these considerations of internal spaces with large curvature, 



33|. 



internal spaces with small curvature have been considered in Ref. 

The position of the brane in the radial direction is not constrained. However, as discussed 
in Sec. (TTJ in the case that the brane is positioned at the center of the disc, only the I = 
KK modes couple to the SM fields. This means that, while the number of kinematically 
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available KK modes is typically of the order of 10 17 for the cases that we consider, in this 
special case this number is effectively reduced to a number of the order of 100. Hence, 
only a negligibly small fraction of the KK modes are effectively available in this case, and 
the signal will be far too weak to be observable at the LHC In order to give a better 
sense for the range of rates possible in the model, it could still be interesting to obtain 
quantitative results for this special case. However, this analysis is complicated by the fact 
that the effective mass splittings between the active KK modes become too large to allow 
us to employ the approximation (|36|) when calculating the total cross section, and thus, 
we have not performed any such calculations. For each fixed set of values for the rest of 
the parameters in the model, we have presented our results for two different values of the 
position of the brane, at r b = r max and at r b = r max /2. 

Note also that, as our calculations are performed to leading order only, we do not take 
final state radiation into account. 



A. p + p—>iet + $ 

For the jet, we have chosen a transverse momentum cut-off p™ m = 750 GeV. The results 
are presented for two different values of the fundamental mass scale, M* = 1.5 TeV and 
M* = 2 TeV. The main background comes from the processes p + p — > jet + Z and p + p — > 
jet + W, with the Z decaying into a neutrino-antineutrino pair and the W decaying into a 



neutrino and a lepton, respectively [34(. In the case of W production, the background can 



be distinguished from the signal if a lepton in the outgoing state is observed. Taking this 



into account, simulations of the background using PYTHIA [35j shows that the background 
from W production is small in comparison to the background from Z production, and hence, 
we have only performed accurate simulations for Z production. 

The differential cross sections dcr/dcos(6 l ) and d<r/dpT,jet are given in Figs. H] and El 
respectively. As a reference, we have also plotted the corresponding cross sections for the 
ADD model with the same value for M*. The cross sections for our model resemble those 
of the ADD model and are of the same order of magnitude. For r b = r max , the result is 
almost indistinguishable from the ADD results for both M* = 1.5 TeV and M* = 2 TeV. As 
expected, the cross sections decrease with increasing M*, while the discrepancy between the 
naive and the truncated cross sections increase with decreasing M*. Also, the signals have 
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the same behavior as the background. For all the demonstrated results, the signals are larger 
than the background, although not by much for M* = 2 TeV. However, for M* = 1.5 TeV, 
the discrepancies between the naive and the truncated cross sections are quite large. 

Note that there is also a difference between results for different values of T5. For both 
values of M* shown, the cross sections are larger for T5 = r max . 

For M* = 1.5 TeV, the integrated cross sections are of the order of 200 fb. Thus, for an 
integrated luminosity at the LHC of 10 fb -1 or 100 fb -1 , the expected number of events is 
of the order of 2000 or 20000, respectively, whereas for M* = 2 TeV, the cross sections are of 
the order of 50 fb, and the corresponding number of events are of the order of 500 or 5000. 




lcos(9)l lcos(9)l 



FIG. 4: The differential cross section for graviton plus jet production with respect to cos(0). The 
solid curves correspond to the parameter value r& = r max and the dashed curves correspond to the 
value Tft = r max /2. For each of these two values, the thick lines are the naive cross sections and 
the corresponding thin lines are the truncated ones. The dotted lines are the corresponding results 
for the ADD model with the same value for M*, and the gray shaded area is the SM background. 
In the left panel, the results are plotted for M* = 1.5 TeV, whereas in the right panel, they are 
plotted for M* = 2 TeV. 

B. p + p^-j + f 

For the photon, we have chosen a transverse momentum cut-off p™ m = 300 GeV. The 
results are presented for M* = 1 TeV and M* = 1.5 TeV. The main background is analogous 
to the background for the jet production process, i. e., coming from the processes p+p — > 7+Z 
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FIG. 5: The differential cross section for graviton plus jet production with respect to pt- The 
solid curves correspond to the parameter value r& = r max and the dashed curves correspond to the 
value Tb = T max /2. For each of these two values, the thick lines are the naive cross sections and 
the corresponding thin lines are the truncated ones. In the left panel, the results are plotted for 
M* = 1.5 TeV and in the right panel, they are plotted for = 2 TeV. 



and p+p — > j+W 



34|. In the same way as in the jet case, the background from Z production 



is dominant, and we have only considered this contribution to the background. 

The differential cross sections da/dcos(6 l ) and dcx/dpr,7 are presented in Figs. M and [3, 
respectively. As for the case of jet production, the cross sections resemble those of the ADD 
model and also have a similar behavior to the background. In this case, the results for 
Tfe = Tmax/2 are very similar to the results for the ADD model. In comparison to the jet 
production case, it is much more difficult to find a region where the effective theory is valid 
and the signal is not much smaller than the background. In particular, in the case that 
M* = 1.5 TeV, the background is much larger than the signal. 

For M* = 1 TeV, the integrated cross sections are of the order of 10 fb. Thus for an 
integrated luminosity at the LHC of 10 fb" 1 or 100 fb -1 , the expected number of events is 
of the order of 100 or 1000, respectively, whereas for M* = 1.5 TeV, the cross sections are 
of the order of 1 fb, and the corresponding number of events are of the order of 10 or 100. 
Comparing to the jet production case, the cross sections are much smaller, about two orders 
of magnitude for M* = 1.5 TeV. 
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FIG. 6: The differential cross section for graviton plus photon production with respect to cos(#). 
The solid curves correspond to the parameter value = r max and the dashed curves correspond 
to the value 75 = r max /2. For each of these two values, the thick lines are the naive cross sections 
and the corresponding thin lines are the truncated ones. In the left panel, the results are plotted 
for M* = 1 TeV, whereas in the right panel, they are plotted for M* = 1.5 TeV. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we have considered a model for large extra dimensions based on hyperbolic 
geometry. More specifically, the internal space has the geometry of a hyperbolic disc with 
constant curvature. This model is in some sense a generalization of the ADD model. Its main 
advantages are that it provides a more satisfactory solution to the hierarchy problem than 
the ADD model, and that constraints from astrophysics on the size of the extra dimensions 
are avoided. 

We have investigated the possible experimental signatures of the model at the LHC. 
The two main reactions that could be relevant for the LHC are p + p — > jet + G and 
p + p — > 7 + G, where G denotes a KK mode of the graviton. As the KK modes are 
extremely weakly interacting, the corresponding amplitudes being suppressed by the Planck 
scale, the gravitons produced in these reactions are not detected, but appear as missing 
energy in detectors. Since there is no known analytic form for the KK spectrum when 
the internal space is hyperbolic, we have employed a combination of approximations and 
numerical investigations to obtain our results. 

We have found that some regions of the parameter space could be probed by the LHC, 
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FIG. 7: The differential cross section for graviton plus photon production with respect to pr- The 
solid curves correspond to the parameter value t& = r max and the dashed curves correspond to the 
value Tb = T max /2. For each of these two values, the thick lines are the naive cross sections and 
the corresponding thin lines are the truncated ones. In the left panel, the results are plotted for 
M* = 1 TeV, whereas in the right panel, they are plotted for M* = 1.5 TeV. 

using the jet production channel. In this reaction, the integrated cross sections are typically 
of the order of 100 fb for the cases that have been studied. For the reaction involving the 
production of a photon, the discovery potential is significantly weaker than in the former 
case. For large M*, the background is much larger than the signal, while for smaller M*, 
the applicability of the effective theory breaks down and the predictions cannot be trusted. 
Also, the integrated cross sections are about two orders of magnitude smaller than for the jet 
production case, typically of the order of 1 fb. For both cases of jet and photon production, 
the cross sections have the same qualitative behavior and are of the same order of magnitude 
as the cross sections for the ADD model with the same value for M*. In fact, the signals of 
our model are in some cases indistinguishable from signals of the ADD model. In addition, 
the parameter space of our model is much larger than that of the ADD model, our model 
having three free parameters with only weak experimental constraints. 

A novel feature of our model, in comparison to the ADD model, is that its physical 
predictions depend on the position of the brane in the internal space. In particular, in the 
case that the brane is placed at the center of the disc, i.e., at r = 0, most of the coupling 
constants between the KK modes of the graviton and SM fields vanish. This has the result 
that in this case the experimental signatures of the model would be far too weak to be 
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observable in collider experiments. 

In conclusion, we have found that the most promising channel for the detection of hyper- 
bolic extra dimensions at the LHC is the production of a KK mode of the graviton together 
with a single hadronic jet. In the case when the fundamental mass scale M* as well as 
the curvature v are of the order of 1 TeV, a solution to the hierarchy problem is obtained. 
Depending on the position of the brane in the radial direction t&, it may also be possible to 
obtain an observable signal at the LHC. 
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